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We discuss the rich scattering features offered by thin and thick plasmonic layers covering dielectric 
nanoparticles and their potential optical applications. The frequency position of scattering dips and peaks 
may be controlled to a large degree using plasmonic layers, which may dramatically vary the total scattering 
signature in the frequency range of interest. We show that complex and exotic spectra may be obtained using 
a single composite nanosphere, including dipole-dipole Fano resonances and electromagnetic induced 
transparency effects. The described phenomena are observable at any angle of observation, in the total 
scattering cross section and for realistic plasmonic materials, since they are based on purely dipolar fields. 
Enhanced field amplitudes are associated with these anomalous scattering features, which may be used to 
efficiently boost weak optical nonlinear effects. Exciting applications of these nanostructures are envisioned, 
such as efficient and tunable sensors, all-optical switches and memories, optical tagging and biomolecular 
imaging. 
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Recently, large interest has been dedicated to the study of plasmonic materials at terahertz, infrared (IR) and 
optical frequencies 1 . At the interface with air or with other dielectrics, these materials support localized 
resonances associated with surface plasmons, of great interest for imaging and energy harvesting applica- 
tions 2 . Specific arrangements of plasmonic and layered nanoparticles may support resonances with sharp spectral 
signatures and large near-field intensities, associated with Fano resonances 3 " 7 and electromagnetic induced 
transparency (EIT) effects 8 " 10 . In general, these extreme resonant phenomena are generated by the strong inter- 
action and coupling between bright and dark resonant modes supported by a complex subwavelength plasmonic 
system. It has been suggested that these exotic scattering responses may have a plethora of practical applications, 
spanning from 'slow light' devices to novel sensing systems. 

In a different context, plasmonic coatings have been used in a nonresonant way to suppress the electric dipole 
scattering from small dielectric or conducting objects, leading to the development of plasmonic cloaking' 
devices 11 " 13 . This concept has been recently experimentally verified at microwave frequencies using plasmonic 
metamaterials 1415 and at IR 16 frequencies. 

It is interesting that plasmonic features can effectively lead to dual scattering phenomena, from very large and 
enhanced scattering signatures 17 " 19 to very low and suppressed response. Hence, exploiting the natural dispersion 
of plasmonic materials we may be able to design a plasmonic nanoparticle configuration supporting low- scatter- 
ing cloaked states followed in close frequency range by highly resonant scattering peaks. In this work we explore 
whether we may combine and control these effects with plasmonic covers, in order to tailor to a large degree the 
complex scattering response of small nanoparticles, spanning Fano, EIT-like, super-resonant and cloaking 
features in a single nanoparticle. 

We limit our analysis to subwavelength spherical nanoparticles, for which the electric dipole response dom- 
inates the overall scattering. Resonant dipolar scattering is generally characterized by a peak with relatively broad 
linewidth. In conventional plasmonic Fano resonances 3 " 6 , this resonance is typically interpreted as the 'bright' 
bonding mode that may be coupled to a narrowband higher-order 'dark' antibonding scattering mode, typically 
an electric quadrupole or higher-order multipole, to generate the typically antisymmetric Fano signature. As we 
have briefly discussed in Ref. 20, however, we may alternatively exploit the 'dark' state associated with the 
plasmonic cloaking scheme 11 to achieve similar sharp scattering responses. This novel mechanism for Fano 
signatures has a different physical nature compared to the conventional antibonding resonant modes, since it 
exhibits lower levels of stored electromagnetic energy and a response that is associated with only dipolar scattered 
fields, and is therefore directly reflected into the total scattering signature of the object. 

In the following, we explore and expand on the concept of Fano -like resonances in realistic composite 
plasmonic nanoparticles, based only on the coupling and interaction of dipolar 'bright' and 'dark' scattering 
states. These dipole-dipole Fano-like resonances have been recently introduced for concentric 2122 , single coated 20 
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and multi-layered 23 geometries. It has been shown that they can 
support strongly confined fields at the nanoparticle's core, which is 
an ideal condition to boost third-order optical nonlinearities to real- 
ize efficient all-optical nanoswitches 20,24 . In the following, single- and 
double -layered nanospheres are theoretically investigated and 
explored. It is shown that they provide unprecedented flexibility in 
the design and manipulation of the total scattering response as a 
function of the thickness and number of plasmonic layers. We fur- 
ther demonstrate that, even with purely symmetric and isotropic 
nanoparticles, it is possible to tailor to a large extent the position 
of zeros and poles of their scattering response and interesting degen- 
erate states may be induced, such as ultrasharp Fano resonances and 
EIT scattering signature, without requiring complex asymmetric 
plasmonic configurations. Some of the concepts presented in this 
report have been briefly discussed in Refs. 20-23 for different pur- 
poses and applied to specific configurations. In Ref. 20 we exploited a 
simple degenerate state supported by a core-shell particle to realize 
all-optical nanoswitches, whereas in Ref. 23 we analyzed the case of a 
four-layer particle supporting a comb-like spectrum with super-scat- 
tering features, which is well suited for optical tagging applications. 
Differently from all these previous works, here we discuss and invest- 
igate more generally the potential of plasmonic cloaks to tailor the 
scattering response and obtain exotic and useful scattering features. 



We also explore their evolution when design parameters and losses 
are considered. These results may have exciting implications for the 
design of advanced photonic and sensing nanodevices. 

Results 

Single layer plasmonic composite nanoparticle. In order to 
understand the possibilities offered by plasmonic cloaks in terms 
of scattering manipulation, we start by considering the basic 
geometry shown in Fig 1(a). It consists of a simple spherical 
nanoparticle composed of a dielectric core (yellow region) with 
relative permittivity 8=2 and radius a coated by a plasmonic shell 
(gray region) of radius a c . In order to gain physical insights into the 
dependence of the scattering properties on the design parameters of 
this nanoparticle, we start by analyzing the ideal lossless scenario in 
which the plasmonic shell is assumed lossless, following a Drude 
permittivity dispersion with parameters 8 c = 8 o0 — fp/\f(f + *Y)]> 
f p =2175 THz, y=0 THz, goo = 5, fitted to the experimental data for 
the real part of silver permittivity 25 , under an e~ iajt time convention. 
Although the lossless case will be useful to gain insights into the 
proposed phenomena, their realistic implementation require 
addressing the issue of losses, which is usually critical in plasmonic 
effects. To assess this important aspect, we also consider in the 
following realistic silver absorption in several cases, corresponding 




frequency f/f p frequency f/f p 

Figure 1 | SCS evolution of a core-shell plasmonic nanoparticle with different aspect ratios, (a) Geometry of a core-shell plasmonic nanosphere. The 
dielectric core has relative permittivity 8=2 and the outer radius is fixed to a c =30 nm in all the considered examples, (b) Contour plot of the normalized 
SCS as a function of the relative permittivity of the plasmonic coating e c /e 0 and the aspect ratio rj c = a/a c . Cloaking (blue lines) and resonant scattering 
(red lines) quasi-static conditions are also shown, (c-f) Normalized SCS versus normalized frequency are shown for lossless (red solid line) and lossy 
(black solid line) plasmonic covers for different aspect ratios: (c) ^ c = 0.95, (d) rj c = 0.7, (e) rj c = 0A and (f) rj c = 0.l. 
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to a collision frequency y=4.35 THz 25 . This will allow assessing the 
sensitivity of the proposed concepts to material losses. 

It is convenient to analyze the scattering response of this spherical 
core-shell geometry by applying the conventional Mie theory 17 . As 
usual, the fields become a discrete sum of orthogonal transverse 
magnetic (TM) and transverse electric (TE) spherical harmonics. 
The n-th order TM scattering coefficient c™ is given as 



where the complex valued quantities U™, V™ can be calculated by 
solving the appropriate determinants 1118 and take real values in the 
lossless case. By duality, one can straightforwardly compute the TE 
complex scattering coefficients cj £ . The total SCS is computed as the 
sum of TM and TE scattering coefficients and it is equal to 

T 2 GO 

scs =iE( 2 " +1 )(l c ™l +V«\\ (2) 

n= 1 

where X 0 is the free-space wavelength. We assume that the composite 
nanoparticle is smaller than the radiation wavelength and has no 
inherent magnetic effects, implying that the TM spherical harmonic 
with n = l dominates Eq. (2), corresponding to electric dipole radi- 
ation. It is usually assumed that, under the condition that only one 
scattering order dominates the total SCS, the scattering spectra may 
not be particularly exotic. However, one may notice that when U™ is 
set to zero, Eq. (1) vanishes and the total SCS may become very small, 
consistent with the cloaking condition described in Ref. 11. Conver- 
sely, when V™ = 0, c™ = — 1 and maximum scattering is obtain- 
ed 18 . If we are able to tailor the position of the zeros of these two 
parameters in the desired frequency range, we may be able to realize 
complex scattering features from such simple symmetric nanoparti- 
cle, even in the quasi-static limit. 

As noticed in Ref. 20, in the long-wavelength regime the quasi- 
static condition to achieve cloaking (U\ =0) depends only on the 
filling ratio rj c = a/a c of the nanoparticle and its material properties: 
rf c (e - Be) (2s c + 1) + (e + 2s c ) (s c - 1 ) = 0, (3) 

Similarly, resonant scattering (V™ = 0) is obtained when the filling 
ratio satisfies 20 

2r\\ (e - Be) (e c - 1) + (e + 2s c ) (Be + 2) = 0 (4) 

We show in Fig. 1 (b) the contour plot of the total SCS of the proposed 
nanoparticle with 8=2 and fixed outer radius a c =30 nm as a func- 
tion of the aspect ratio r\ c and of the relative permittivity s c of the 
plasmonic coating (brighter colors correspond to larger SCS). In this 
panel, we consider the lossless scenario only. The dispersive permit- 
tivity e c increases monotonically with frequency, as dictated by 
Kramers -Kronig relations and indicated by the horizontal dashed 
arrows in Fig. 1(b). The dispersion of the solutions of Eqs. (3)-(4) 
is also plotted in Fig. 1(b), as the solid blue (cloaking) and red (res- 
onant scattering) lines, respectively. 

In Figs. l(c)-(f), we present the evolution of the normalized total 
SCS spectra (solid red lines) for different aspect ratios, as indicated by 
the different horizontal dashed lines in Fig. 1(b). The total SCS of a 
bare dielectric object (dashed blue line) with same outer radius as the 
considered nanoparticle is also plotted in Figs. 1(c) -(f) for compar- 
ison. It is interesting to see how the intersections of the dashed lines 
in Fig. 1(b) with the blue and red dispersion curves (zeros of U™ and 
V™, respectively) are directly reflected in dips and peaks in the total 
SCS. When the location of these zeros gets closer in frequency, 
anomalous scattering features are observed. Figure 1(c) shows the 
scattering response of a nanoparticle with large aspect ratio 
(f/ c =0.95), which corresponds to a thin plasmonic shell. In this case, 
a rather broadband low scattering (cloaked) region is observed at 
low frequencies, similar to the conventional plasmonic cloaking 



scheme 11 . At higher frequencies, a second cloaking condition lies 
in close proximity to a resonant state [see Fig. 1(b)], producing a 
very peculiar asymmetric response. This is a first example of degen- 
erate state between cloaking and resonant scattering, which resem- 
bles the shape of a Fano resonance. We stress that this is rather 
different from conventional Fano resonances, in which different 
orthogonal scattering orders interfere within the angular spectrum. 
Here, the same scattering order (electric dipolar) is responsible for 
this anomalous scattering response. Figure 1(d) shows the SCS spec- 
trum for a plasmonic nanoparticle with aspect ratio rj c =0.7. In this 
case, two distinct frequency regions of suppressed scattering are 
obtained and the nanoparticle can be effectively transparent in two 
different frequency bands 12 . However, a high resonant scattering 
peak always exists between the two scattering dips. This is a general 
constraint when considering causality and passivity limitations of 
plasmonic materials, which prevents the possibility of realizing a 
passive broadband cloaking device based on merging two cloaked 
states along the frequency axis 12,23 . Always a scattering peak is 
expected to arise between the two dips, consistent with Foster's 
reactance rule for the response of passive circuit networks 26 ' 27 . 

For lower aspect ratios the scattering signature of the proposed 
plasmonic nanoparticle becomes even more intriguing and complex. 
The SCS for rj c =0A is shown in Fig. 1(e) and again a relative broad 
cloaking frequency region is observed, now in the higher frequency 
range. More interestingly, an additional asymmetric scattering res- 
onance is obtained at lower frequencies, due to the interference 
between the closely spaced cloaking dip and scattering peak. As in 
panel (c), the merging of the cloaking and resonant states gives rise to 
a dipole-dipole Fano-like resonance 20 " 23 , which can have extremely 
narrow linewidth for even lower aspect ratios, as observed in Fig. 1(f) 
for rj c =0.1. The linewidth is ultimately limited by losses in the mate- 
rials and quantum effects that bound how small we can make the 
inner core without affecting the resonant state. To quantitatively 
assess the sensitivity to realistic material absorption, we also show 
the SCS of the composite particle when the real losses are considered 
[black curves in panels (c)-(f)]. It is evident that absorption affects 
the scattering signature, by reducing the scattering excursion, but the 
general behavior and the most interesting features are maintained. It 
is interesting that the asymmetric resonances can retain a very sharp 
linewidth, except when the cloaking and resonance states are too 
close in frequency, as in panel (f). This remarkable robustness to 
losses is due to the purely dipolar nature of the proposed effect. 

It is important to stress that, different from conventional plasmo- 
nic Fano-resonances 3 based on the interaction of 'bright' dipolar with 
'dark' higher-order modes, this dipole-dipole Fano-like response can 
be observed for any angle of observation and it is detectable in the 
total SCS signature, due to the isotropy of the plasmonic nanoparticle 
and the dipolar nature of both bonding and antibonding states 
involved. 

To further study these interesting phenomena, we explicitly cal- 
culate the field distributions for rj c = 1/3 at specific frequencies of 
interest for the lossless geometry. The normalized SCS versus wave- 
length is plotted in Fig. 2(b) in the lossless case. The points I and IV in 
Fig. 2(b) are the usual, moderately broad cloaking and resonant 
scattering states of the plasmonic nanoparticle, respectively. 
Instead, the points II and III in Fig. 2(b) correspond to an extremely 
narrowband asymmetric resonance with an almost 50 dB abrupt 
change between cloaking and high scattering states. As explained 
in the previous paragraphs, the asymmetric Fano-like scattering res- 
ponse is a signature of the degenerate state between these two distinct 
dipolar phenomena. 

The corresponding electric field distributions on the E plane are 
plotted in Figs. 2(c) (snapshot in time of the E x component, parallel 
to the incident electric field) and (d) (magnitude of the total electric 
field) for all these frequency points. In addition, the arrows in panel 
(d) represent the power flow (time-average Poynting vector) 
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Figure 2 | Field plots of a core-shell plasmonic nanoparticle. (a) Geometry of the plasmonic nanosphere, with outer radius a c =30 nm, aspect ratio 
r\ c = 1/3 and 8=2. b) Normalized SCS versus wavelength in the lossless case, (c) Corresponding electric field plots for the real part of the E x component 
(snapshot in time) and (d) the amplitude of the total electric field, in the E plane at the broad cloaking dip (I), the narrow Fano-like resonant peak (II), the 
narrow Fano-like cloaking dip (III) and the broad resonant peak (IV). Blue arrows in (d) represent the corresponding power flow (time-average Poynting 
vector) distributions. The impinging plane wave travels from bottom to top with electric field linearly polarized along x. Amplitude levels in all panels are 
normalized to the impinging wave. 



distribution around the nanoparticle. The overall field distributions 
at points II and IV show a clearly dipolar resonant response. The 
corresponding power flow around the particle is very complex, since 
the SCS is much larger than the physical cross section of the particle 
and typical vortices and saddle points are observed in the vicinity of 
the particle 28 ' 29 . The fields in the core are very much enhanced at 
these two wavelengths, as expected for a resonant state. At the other 
two wavelengths of interest, points I and III, a dual scattering res- 
ponse is observed: very low scattering, i.e., cloaking, is achieved right 
around the nanoparticle and the impinging plane wave fields are 
effectively restored, as if the particle was transparent. By comparing 
the two cloaked cases, however, a crucial difference may be observed: 
at the broad cloaking frequency (point I), the field amplitude is also 
very low inside the core, comparable to the impinging levels; con- 
versely the fields are strongly enhanced and confined at the particle's 
core at the Fano dip (point III). This is due to the close interaction of 
the cloaked and resonant states in this latter case, which produces the 
exotic situation in which the fields are almost unperturbed outside 
the cloak, but strongly enhanced inside it. This phenomenon can lead 
to efficient sensors and enhanced nonlinear effects 20 . It is remarkable 
that in a small frequency window (points II and III), the field distri- 
bution outside the nanoparticle can be dramatically altered from 
invisible to a bright scattering state, while the interior fields are 
drastically enhanced. These interesting properties, combined with 
enhanced third- order nonlinear effects, may potentially lead to 
highly efficient all-optical nanoswitches and memories 20 . 

Double-layer plasmonic composite nanoparticle. After having 
established that a simple plasmonic layer can produce exciting 



scattering features, we show in this section that even more 
complicated scattering signatures may be designed and tailored 
using multilayered cloaks for a single, isotropic subwavelength 
nanoparticle. As we show in the following, multiple Fano-like 
resonances and other complex scattering features may be produced 
using two different plasmonic coating layers, further enriching the 
nanoparticle's scattering spectrum. The geometry of this configur- 
ation is shown in Fig. 3(a): the dielectric core has a relative permit- 
tivity 8i=2 and the outer radius of the composite nanoparticle is 
fixed at a c =50nm. The dispersion of the first plasmonic layer 
follows the same Drude permittivity, and we consider both the 
lossless (y=0 THz) and lossy scenarios (y=4.35 THz). The second 
plasmonic layer has a relative permittivity s 3 with similar Drude 
dispersion, but a plasma frequency decreased by A/p = 325THz. 
This may be achieved, for instance, by modifying the doping level 
of some plasmonic semiconductors, similar to the case of four-layer 
plasmonic nanoparticles based on aluminum- doped zinc oxide 
(AZO) semiconductors 30 ' 31 proposed in Ref. 23. 

In this case, the quasi- static condition for resonant scattering 
(V™ = 0) may be calculated using Mie theory: 



2f]\ (ei - e 2 ) [(e 2 - e 3 ) (e 3 +2)+n\ (2s 2 + e 3 ) (e 3 - 1)] + 
r}\ ( fil + 2e 2 ) [2f]\ (e 2 - e 3 ) (e 3 - 1) + (£2 + 2s 3 ) (e 3 + 2)] = 0, 



(5) 



where rj l =(a + ti n )/(a + t in + t out ) is the aspect ratio of the inner 
plasmonic layer, with radius a+t in , over the outer plasmonic layer, 
with radius a + t in + t out =a c . Furthermore, rj 2 = a /{a + t in + t out ) is the 
aspect ratio of the dielectric core with radius a over the outer 
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Figure 3 | SCS of a double-layer plasmonic nanoparticle with different inner layer thicknesses, (a) Geometry of the double-layered plasmonic 
nanosphere. The dielectric core has relative permittivity Si=2 and the outer radius is a c =50 nm. The outer layer thickness is fixed to t out =5 nm. (b), (c) 
and (d) show the normalized SCS for lossless (red solid line) and lossy (black solid line) plasmonic covers for three different inner layer thicknesses: 
ti n =l nm, tf n =10 nm and t{ n =20 nm, respectively. 



plasmonic layer. The parameters t in and t out are the thicknesses of the 
inner and outer plasmonic shells, respectively. In a similar fashion, 
the cloaking condition (17™ = 0) can be evaluated in the long- wave- 
length limit as 

l/l(ei-B2)[2(B2-fi3)(fi3-l)+l??(2fi2 + fi3)(2fi 3 + l)] + 

, r , , (6) 

rj\ (ei + 2s 2 ) [n\ (e 2 - s 3 )(2s 3 + 1) + (e 2 + 2s 3 ) (e 3 - 1)] = 0. 

As in the previous section, these formulas can accurately describe the 
cloaking and resonant scattering behavior of the double-layer par- 
ticle in the quasi-static limit of interest here. 

In Fig. 3, we show the scattering response of this double-layer 
structure varying the thickness of the inner layer. While this thick- 
ness is modified, the outer radius of the nanoparticle is kept constant, 
i.e., we shrink the dielectric core, for fair comparison. The outer layer 
thickness in this study is fixed at t out =5 nm. In this case, three 
distinct scattering peaks and dips may be induced in the scattering 
signature, visible in Figs. 3(b), (c) and (d), where the normalized SCS 
is plotted versus impinging wavelength for three different inner layer 
thicknesses: ^=1 nm, t in = 10 nm and t in =20 nm, respectively. 
Figure 3(b) clearly demonstrates two consecutive narrowband 
Fano-like resonances indicated by the labels CI, Rl and C2, R2, 
respectively. Each Fano-like feature is a degenerate state between 
neighboring cloaking and resonant scattering, generalizing the 
results of the previous section. We notice in the following panels 
that, as the inner layer thickness is increased, the scattering spectrum 
evolves to an interesting "inverse EIT" signature around 2=360 nm 
(Panel (c) in Fig. 3), due to the strong interaction and coupling 
between the narrow linewidths of the two distinct Fano-like reso- 
nances. This can be considered a triple degenerate state between a 
resonant mode and two cloaking states. These effects are mostly 
preserved after adding realistic level of losses in the plasmonic layers 
[black curves in Fig. 3]. Finally, in Fig. 3(d), we can see how the 
inverse EIT profile collapses to an ultrasharp Fano resonance, as 
the thickness of the inner layer is further increased. However, in this 
case absorption strongly affects this ultranarrow resonant feature. It 
is exciting to observe the flexibility with which we can engineer the 
scattering spectrum of a simple symmetric composite nanoparticle 
by simply moving the zeros of U™ and V™ by varying the geometry 
according to Eqs. (5) -(6). 



So far, we have varied only the inner layer thickness. Analogously, 
the thickness of the outer layer may also be used as a tool to manip- 
ulate the scattering response, as we consider in Fig. 4. The figure is 
analogous to Fig. 3, but in this case the inner layer thickness is fixed to 
t in =5 nm and the scattering response is analyzed as the thickness of 
the outer layer is increased. As seen in panels (b), (c) and (d) of Fig. 4, 
also in this case the two Fano-like resonances are modified in a 
peculiar way as the geometry of the nanoparticle is changed. In 
particular, it is interesting that in Fig. 4(c) the two Fano resonances 
"unfold" and the scattering signature assumes a symmetric EIT-like 
shape around 350 nm, unaffected by the realistic plasmonic losses, 
with interesting potential for multiband sensing. 

It is clear from these results that dual-layer cloaks offer exciting 
possibilities to manipulate the scattering spectrum of small nanopar- 
ticles, even if limited to just the dominant, purely dipolar scattering 
order. By varying the position of the zeros of U™ and V™, we are 
able to position cloaked and resonant states across the scattering 
signature spectrum and tailor the overall response of the nanoparti- 
cle to a large degree. This is in some sense similar to the way filter 
responses are designed in circuit theory by manipulating the position 
of zeros and poles of their transfer function 26,32 . These results are 
summarized in Fig. 5, which shows the frequency variation of the 
different cloaking and resonant scattering conditions as a function of 
the inner and outer layer thicknesses. Comparing panel (a) (inner 
layer variation) and (b) (outer layer variation) of Fig. 5, it is inter- 
esting to note that the Fano-like resonance composed by the points 
C2 and R2 is more sensitive to the change of the outer layer thickness 
than to the inner layer one. Conversely, the Fano-like resonance 
arising at points CI and Rl is affected by the thickness of the inner 
layer. The branches C3 and R3 are influenced by a change of both 
inner and outer layer thicknesses, consistent with the results in Figs. 3 
and 4. Therefore, each cloaking- resonance pair is sensitive in a dif- 
ferent way to the geometrical parameters of the double-layer nano- 
particle and this can be used to engineer the spectrum in a variety of 
ways. Figure 5 also provides an insight on how the Fano-like response 
is physically established: it is evident that the cloaking-resonance pair 
Cl-Rl is due to the interaction of two surface plasmons on the two 
sides of the inner plasmonic layer, whereas the pair C2-R2 is deter- 
mined by the plasmons on the outer plasmonic layers. We have 
verified by inspecting the field distributions at the two wavelengths 
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Figure 4 | SCS of a double-layer plasmonic nanoparticle with different outer layer thicknesses, (a) Geometry of the double-layered plasmonic 
nanosphere, similar to Fig. 3. Here the inner layer thickness is fixed to f jn = 5 nm. (b), (c) and (d) show the normalized SCS for lossless (red solid line) and 
lossy (black solid line) plasmonic covers for three different outer layer thicknesses: t out = 1 nm, t out = 10 nm and t out =20 nm, respectively. 



that this is indeed the case (not reported here for brevity). Finally, the 
C3-R3 pair, sensitive to both inner and outer layer thicknesses, is 
based on less confined plasmon modes, which is reflected in less 
sharp spectral features. It is evident that the addition of a second 
plasmonic layer provides more degrees of freedom to the design and 
tailoring of the scattering response. As demonstrated here, we may 
achieve significantly complex scattering signatures, such as multiple 
Fano-like resonances which can also evolve in EIT-like responses. 

Discussion 

In this Report, we have demonstrated an unprecedented capability to 
manipulate and tailor the scattering response of a composite nano- 
sphere with subwavelength size by using multilayered plasmonic 
covers. We have first investigated the simple core -shell geometry 
and we have shown interesting degenerate states between cloaking 
and resonant scattering, which give rise to ultrasharp dipole-dipole 
Fano resonances in the scattering spectrum. The addition of a further 
plasmonic layer leads to even more complex scattering responses, 
beyond conventional Fano resonances, such as double and triple 
degenerate states, resembling an EIT-like response. Since the inter- 
ference between the different states is very sensitive to the geometry, 
a plethora of different shades between Fano and EIT lineshapes can 



be obtained by simply modifying the aspect ratios of the multi- 
layered nanoparticle. More layers may directly lead to even higher 
flexibility in terms of scattering manipulation, allowing exotic scat- 
tering signatures, such as the "Fano comb" spectrum discussed in 
Ref. 23. As expected, when material losses are considered, the Fano- 
like resonance linewidth slightly broadens and the difference 
between scattering peaks and dips becomes less pronounced. 
However, our results show that the dipolar nature of the involved 
phenomena leads to high robustness to realistic level of plasmonic 
material losses and, as a result, the peculiar shape of the scattering 
spectrum is preserved even after absorption is taken into account, 
which is an important feature towards the practical implementation 
and experimental verification of the proposed concepts. 

These findings may have interesting applications in several fields 
of optics, ranging from signal processing to optical tagging, biological 
sensing and imaging. The proposed multilayered plasmonic nano- 
particles may realize optical narrowband bandpass nanofilters 33 ' 34 
whose response may be tailored with great flexibility. Moreover, 
since some of the scattering features outlined here, particularly those 
associated with surface plasmons supported at the outer interface, 
may be very sensitive to the background medium, these results may 
open interesting scenarios for biological and chemical sensing 6 . 




0 5 10 15 20 25 0 5 10 15 20 

Inner Layer Thickness (nm) Outer Layer Thickness (nm) 



Figure 5 | Frequency variation of cloaking and resonant scattering states versus the thickness of the plasmonic layers. Variation in wavelength of the 
cloaking dips (CI, C2, C3 blue lines) and resonant peaks (Rl, R2, R3 red lines) as a function of the thickness of the nanop article's (a) inner layer and (b) 
outer layer. In both panels the geometry is shown in the inset and the layer under investigation is highlighted as a red striped surface. The vertical dashed 
lines correspond to the cases reported in panels (b), (c), (d) of Figs. 3 and 4, as indicated by the labels. 
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Conversely, other scattering states, localized in the internal interfaces 
of the nanoparticle, may be very robust to the background medium 
and, for this reason, may be ideal to develop optical tagging plat- 
forms 23 independent of the background, with several potential appli- 
cations. For example, the extremely sharp narrowband scattering 
response, discussed in some of the previous examples, is very well- 
suited for multicolor labeling in optical bio -imaging techniques 35,36 . 
The possibility to cover the visible and near-infrared frequency range 
with a large number of narrow resonance peaks may allow detecting 
multiple molecular targets at the same time, with great benefits, e.g., 
for tissue diagnosis. The flexibility of our approach for the manip- 
ulation of the scattering signature of subwavelength nanoparticles 
appears ideal to meet the requirements of all these applications. 

Methods 

All the results presented in this Report are obtained using analytical calculations 
based on full-wave Mie theory for scattering of spherical objects (for further details 
and physical insight, see, e.g., Ref. 17). As usually done, the impinging plane wave and 
the scattered fields are expanded in TE and TM spherical harmonics in a spherical 
coordinate system centered with the nanoparticle. By applying proper electromag- 
netic boundary conditions on each spherical interface, the scattering coefficients for 
TE and TM waves are obtained, which uniquely determine the fields everywhere. 
Total scattering cross section and field distributions are computed using Eq. (2) and 
related formulas, truncating the summations after reaching proper numerical 
convergence. 
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